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Isomerization between the cis and trans conformers of the Si state of acetylene is studied using a reduced 
dimension DVR calculation. Existing DVR techniques are combined with a high accuracy potential energy 
surface and a kinetic energy operator derived from FG theory to yield an effective but simple Hamiltonian 
for treating large amplitude motions. The spectroscopic signatures of the Si isomerization are discussed, with 
emphasis on the vibrational aspects. The presence of a low barrier to isomerization causes distortion of the 
trans vibrational level structure and the appearance of nominally electronically forbidden A ^ Ai -s— X 
transitions to vibrational levels of the cis conformer. Both of these effects are modeled in agreement with 
experimental results, and the underlying mechanisms of tunneling and state mixing are elucidated by use of 
the calculated vibrational wavefunctions. 
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I. INTRODUCTION 

The history of the A state of acetylene is full of spec- 
troscopic surprises. Beginning with the discovery that 
acetylene changes shape from linear to trans-bent upon 
electronic excitation^r— , the A -i^ X spectrum has exhib- 
ited a string of unexpected phenomena, including axis 
switching^, triplet perturbations-"^, as well as unusually 
strong Coriolis interactions and Darling-Dennison reso- 
nance between the bending vibrations^. Most recently, 
the comprehensive assignment of the low energy vibra- 
tional structure has led to the identification of several 
"extra" levels. These appear in the spectrum near 3000 
cm~^ above the ground vibrational state of the trans con- 
former, just as the effective Hamiltonian (H°*^) models 
developed at lower energy begin to break down^*^. 

The most probable explanation for these "extra" lev- 
els is that they belong to the cis conformer of the Si 
stateii. Ah initio calculations have long predicted a cis 
minimunii^"— , but experimental confirmation has been 
hard to come by. This is because the Si-Sq electronic 
transition is forbidden in Ci^u geometries, and any spec- 
troscopic observation of the cis conformer is therefore ex- 
pected to manifest itself in the form of weak transitions or 
indirect effects, necessarily resulting from vibronic mech- 
anisms. 

The purpose of this paper is to investigate theoreti- 
cally the spectroscopic consequences of the low barrier 
isomerization process on the A state potential surface. 
In general, how will the presence of a second distinct 
accessible minimum distort the vibrational energy level 
patterns? One detail of particular interest is whether Si 
cis states could appear at the observed energies with the 
intensities observed for the "extra" levels. The calcula- 
tions in the literature characterize the relevant station- 



ary points on the potential energy surface (PES), but 
cannot address these questions because of the delocal- 
ized nature of the wavefunctions involved in large am- 
plitude motion and the importance of effects like tun- 
neling. Large amplitude dynamics on a reduced dimen- 
sion Si surface have been studied previously by spectral 
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, but those calculations focused mainly 



on the well-known Franck-Condon active progressions in 
the CC stretching mode, vi, and the frans-bending mode, 
W3, while revealing little about the weak or forbidden 
bands that encode the isomerization dynamics. 

The method we choose for our treatment uses a 
reduced dimension Discrete Variable Representational. 
Since we are investigating an isomerization process, the 
calculation must be able to treat large amplitude mo- 
tions that span the two geometries. DVR methods are 
well-suited to this type of problem because the basis func- 
tions are not tied to a single center as they usually are in 
a variational calculation^®. Since the half- linear transi- 
tion state^"^ and the Franck-Condon active vibrations are 
all planar, we can perform a reduced dimension calcula- 
tion and still expect to capture the essential features of 
the experimental spectrum. As described in detail below, 
the calculation is performed in the three dimensions that 
encompass the CC stretch as well as the trans and cis 
bending vibrations. 

An important ingredient in our calculation is a high 
accuracy potential energy surface for an excited elec- 
tronic state. This work represents the first application 
of EOM-CCSDTlS methods to the PES of a polyatomic 
molecule, and we believe there is great promise in apply- 
ing such methods to the study of electronically excited 
states. In fact, we find that this level of theory is neces- 
sary to obtain quantitatively acceptable agreement with 
experiment for this system. 
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We obtain from our DVR calculation the vibrational 
eigenstates of the full Si potential surface in the three di- 
mensional coordinate space. Since the basis functions are 
not defined relative to one specific geometry, the eigen- 
functions are not predetermined to belong to either the 
cis or trans conformer, although they naturally divide 
themselves in this way at low energy. We therefore ex- 
pect that the calculated results will contain all of the 
effects we desire to model: the possible existence of cis 
levels interspersed among those of the trans conformer, 
and any mixings between them that arise from tunnel- 
ing through the isomerization barrier. Other signatures 
of the isomerization, such as the distortion of the trans 
level structure, should also be satisfactorily reproduced. 

Finally, we note that this method could be gener- 
ally applicable to other isomerizations, even in larger 
molecules. It allows for a selective treatment of only the 
few coordinates relevant to the minimum energy isomer- 
ization path, but is still based on a simple and easily 
understandable Hamiltonian. Its foundation on estab- 
lished DVR methods provides computational efficiency 
and accuracy. 



II. METHODS 

In this section we describe the elements of our 
method for reduced dimension vibrational DVR calcu- 
lations. Ours is not the first such calculation, even on 
acetylene22i2i, but our approach differs from those of pre- 
vious authors. We specify significant departures from 
earlier reduced dimension DVR calculations on acetylene 
where appropriate. 

An important consideration in a reduced dimension 
calculation is that the qualitative and quantitative ac- 
curacy of the results are frequently balanced against the 
effort required to perform the calculation. In this work 
we have gone to some lengths in order to obtain quantita- 
tively accurate results, but very useful qualitative results 
can be had much more easily. Even for spectroscopists, 
the pattern and assignments of the levels can be valu- 
able even if quantitative agreement with experiment is 
poor. We demonstrate this point later on in connection 
with the quality of the potential energy surface and the 
number of dimensions included in the calculation, which 
are the two factors that have the greatest impact on the 
computational cost of the calculation. 

It is appropriate to begin our discussion of the details 
by mentioning that we desire a matrix representation of 
H = T + V, where T and V are expressed in the same 
basis. We will examine the two parts of the Hamiltonian 
separately, but first we must digress briefly on the subject 
of coordinate systems. 



A. Molecule-Fixed Coordinate Systems 

The choice of coordinate system for a calculation in 
the molecular frame is a complicated one. The problem 
of separating vibration and rotation is exacerbated when 
large amplitude motions are possible^-, and even remain- 
ing in the center of mass frame may not be trivia P^i^^ for 
molecules larger than triatomics. For acetylene, several 
different coordinate systems have been used in the lit- 
erature for various types of calculations^^i^. We find 
the internal coordinates of FG theory^^ to be extremely 
convenient, and recommend them generally. The poten- 
tial energy surface is compactly represented in these co- 
ordinates, especially at low energy, and the G matrix 
elements^i^l derived from them provide a simple way to 
represent the kinetic energy operator. 

In reduced dimension calculations using internal coor- 
dinates, one departure that may be necessary from the 
canonical set of such coordinates for a given molecule 
concerns their domains. Since the usual set of 3A^ — 6 
coordinates is chosen to uniquely specify all possible ge- 
ometries (leaving aside linearity for the moment), it may 
be necessary to change the domain of one or more of the 
angles. For example, since we do not include the usual 
torsional coordinate, r, in our planar acetylene calcula- 
tion, it is necessary to increase the domain maximum of 
at least one of the ZCCH valence angle bends from tt to 
27r in order to allow both cis and trans geometries. (In 
this paper r is defined as zero, except in Table |T] where 

T = TT.) 

B. The Potential Energy Surface 

The hallmark of all DVR calculations is the simplicity 
of the potential energy matrix elements. More gener- 
ally, any function of coordinates is a diagonal matrix, 
{di\f{x)\9j) « f{xi)5ij, where the {xi} are the pseu- 
dospectral grid points of the DVR basia^. Once a re- 
duced set of coordinates has been chosen, only two things 
remain to be done to express V in the DVR basis: find- 
ing the appropriate level of ab initio theory at which to 
calculate the potential energy surface, and deciding how 
to treat the discarded degrees of freedom when the active 
coordinates are varied. 

Previous theoretical investigations of the Si surface 
encountered difficulties with the harmonic frequencies 
of the trans structure, especially the lowest frequency 
modes, U4 and vu^^. We report here EOM-CCSDT 
harmonic frequencies and compare them to EOM-CCSD 
values and experimental data in TableUl Since we exclude 
V4, the CCSD frequencies seem satisfactory, but still the 
improved agreement at CCSDT is obvious. Indeed, we 
were able to obtain qualitatively useful results with a 
CCSD potential surface, but quantitative agreement re- 
quires CCSDT. In light of this, a reasonable full set of 
harmonic frequencies appears to be a good criterion for 
selecting a level of ab initio theory. The relatively poor 
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TABLE I. Geometries and harmonic frequencies for A ^ 
^^C2H2: Comparison of ah initio methods 







CCSDT 


Expt.^^-^ 




3107.7 


3053.5 


(3004 


102 


1471.4 


1424.6 


1410.26 




1106.6 


1106.1 


1070.34 


UJi 


614.6 


761.5 


764.7 


W5 


3083.7 


3031.0 


(2914^^ 


0^6 


745.8 


776.5 


772.5 


ZCCH 


123.64 


122.09 


122.48 


Rcc 


1.3575 


1.3743 


1.375 


Ren 


1.0907 


1.0963 


1.097 


To 


43830 


42334 


42197.57 



Frequencies in cm ^, angles in degrees, and bond lengths in A. 
These harmonic frequencies are less well determined because 
their overtones have not yet been observediffi 

quantitative agreement with experiment obtained in a 
full dimensional calculation of the vibrational fundamen- 
tals using a CCSD surface^^ corroborates this conclusion, 
as our reduced dimension CCSD vibrational fundamen- 
tals for V3 and ve are very similar to their results. 
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FIG. 1. A slice of the Si PES with -Rcc ^^nd -RcH at their trans 
equilibrium values. The cis minima are located in the lower left 
and upper right, and the trans minima are in the upper left and 
lower right. The energy scale is in cm~^ and the contour spacing 
is lO-^ cm~l. 

The simplest approach for dealing with the discarded 
degrees of freedom is to fix them at their equilibrium 
values. While it is impossible to make absolute gen- 
eralizations due to the heterogeneous decomposition of 
the normal modes into the internal coordinates, our ex- 
perience is that this approximation leads to computed 
vibrational fundamentals within at worst 10% of the ex- 
perimental values, but usually significantly better. We 
would have been satisfied with this approximation except 
for the Fermi resonance that manifested itself between 



2ragm g^jjj 2"-i6'""'"^ in early calculated results, which 
precluded comparisons with effective Hamiltonian fits to 
experimental data. We therefore replace i?cc with the 
normal coordinate for V2 from the harmonic frequency 
calculation, which provides a prescription for varying the 
CH distances as the CC bond length changesi^ We did 
also consider vibrational adiabatic potentials, which are 
commonly used in the literature^Si^i. However, given the 
limited return and the enormous computational cost of 
minimizing the energy at every point on the entire three 
dimensional grid, we elected not to pursue this option. A 
fitted surface might be more suitable for addressing this 
particular issue. 

C. The Kinetic Energy Operator 

The kinetic energy operator presents several difficul- 
ties. The form of the Laplacian in the chosen coor- 
dinate system must be determined, which can be very 
tricky22i2ii^j especially since singularities in the kinetic 
energy are possible. Furthermore, expressing the Lapla- 
cian in a DVR basis can lead to a non-Hermitian H. 
The kinetic energy operator matrix elements are more 
complicated than those of V, but we merely need to 
add matrix representations of differential operators to the 
trivial matrix representations of functions of coordinates 
already in our repertoire. Depending on the underly- 
ing basis, DVR matrix elements of differential operators 
can be calculated straightforwardly either via basis trans- 
formations to and from a finite basis representatioiii^, 
e.g. a harmonic oscillator, or by explicit formulas^, most 
commonly based on Fourier functions^^"— . The selec- 
tion of an appropriate DVR basis consists primarily of 
choosing one whose underlying basis functions have the 
same boundary conditions as the eigenfunctions to be 
obtained'^''. 

We use the G matrix elements to obtain the kinetic 
energy operator, taking care to follow Podolsk}*^, and 
also using an explicitly Hermitian forir4S of each term, 
such that 

where g — det |G| and the Si are the internal coordi- 
nates. It is important to note that although the tradi- 
tional FG matrix solution invokes the approximation of 
infinitesimal displacements, the coordinates themselves 
and their G matrix elements are valid for motions of any 
amplitude'^''' . We simply represent the functions of coor- 
dinates in the Gij as diagonal matrices instead of evalu- 
ating them at a specific geometry. In our planar acetylene 
calculation we take care to use the kinetic energy matrix 
elements for the bending of a linear molecule^!, since the 
torsional angle is not well behaved when planar cis-trans 
isomerization can occur. 

Finally, for consistency with our use of the V2 normal 
coordinate when calculating the PES, we also redefine 
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Rcc as the CC stretch component of the V2 normal co- 
ordmate, according to Qk = X]* {^~^)kt '^^^'^^ Q' 
and L are the mass- weighted normal coordinate, an inter- 
nal coordinate, and the transformation that diagonalizes 
GF, respectivel y^^ i'*^ . The Rch contributions are still 
neglected. 

D. Specifics and Practical Considerations 

The specifics of our treatment of Si acetylene are 
as follows. We take a subset of the internal coordi- 
nates, consisting of the two ZCCH valence angle bends 
(0 < {01,02} < Stt) and the C-C bond stretch, which 
correspond loosely to {v2, v^^vq}. Our DVR basis is then 
the direct product basis 01 (8) 02 (8) Rcc , and we use the 
sine DVR basis'^'' for the angles, but a harmonic oscillator 
basis for the bond stretch. The final results were gener- 
ated using 120 grid points for each of the two bend bases, 
and 43 points for the stretch basis. The pseudospectral 
grid points and differential operator matrix elements are 
available from analytical formulas for the sine basis, and, 
for the harmonic oscillator basis, from the diagonaliza- 
tion of the coordinate matrix and the application of the 
resulting transformation matrix to other operators ex- 
pressed in the harmonic oscillator basis. We then con- 
struct T according to Eq. [TJ We calculate the potential 
energy for 10 geometries along each coordinate at regular 
intervals spanning {60° < 0i < 300°, 60° < 02 < 180°, 1 
A < Rcc < 2 A}. Ab initio calculations were performed 
with the CFOUR program system^i^, using the EOM- 
CCSDT method and the NASA Ames ANOl basis set. 
The elements of V are found by interpolating the po- 
tential energy surface at the DVR grid points, after dis- 
carding grid points that lie outside the original domain 
of the ab initio surface. This reduced the size of the basis 
mentioned above to 80x80x34. 

In practice, the construction and diagonalization of H 
are accomplished by a basis set contraction^^. The DVR 
H for the 2D 0i(8)02 space at every value of Rcc is diago- 
nalized, and the eigenvectors above a chosen cutoff energy 
are discarded. (For the results presented here this cutoff 
energy was 10,000 cm~^, with a minimum of 35 vectors 
retained per 2D block.) The resulting rectangular trans- 
formation matrices are used to compress the remaining 
blocks of the 3D H, which are off-diagonal in i?cCi prior 
to the final diagonalization. Applying the rectangular 
transformation matrices in reverse transforms the final 
eigenvectors back to the grid point representation. 

III. RESULTS 

Using the method described here, vibrational eigen- 
states of Si acetylene could be obtained up to energies 
exceeding 15,000 cm~^. For the most part we will limit 
the discussion to states up to about 5,000 cm~^. In gen- 
eral we will forego itemized examination of the calcu- 



TABLE II. DVR vibrational fundamentals for A M„ and 
A '^CaHa 





A 




A ' 


A2 




Calc. 


Expt.^ 


Calc. 


Expt.-ii 




1414.65 


1386.9 


1489.61 






1033.6 


1047.55 


789.56 






780.35 


768.26 


588.35 


(56511 



Frequencies in cm . 
^ Estimated from the energy difference between two observed 
combination bands. 



TABLE III. Ab initio geometry and harmonic frequencies for 

A '^C2H2 





2997.14 


(JJ2 


1583.22 


UJ-i 


806.09 


UJ4 


817.5 


LJ5 


2941.81 




571.62 


ZCCH 


132.62 


Rcc 


1.3423 


Rch 


1.0983 


To 


45155 



Frequencies in cm ^, angle in degrees, and bond lengths in A. 

lated states and instead focus on more broad agreement 
with experiment and predicted trends, in keeping with 
the goals stated in the beginning of Section |lTl Never- 
theless, it is worth mentioning that the calculated vibra- 
tional fundamentals are within 2% of the experimental 
values (Table HI]), despite the neglect of the other dimen- 
sions. (The results of an ab initio harmonic frequency 
calculation for the cis geometry are given for reference 
in Table Uni) 

A. Symmetry of the Reduced Dimension Eigenstates 

Before entering into a more detailed discussion of the 
calculated states, it is necessary to work out the connec- 
tions between the symmetries that exist in the coordinate 
space of the calculation and the true symmetries of the 
molecule, so as to be able to interpret the wavefunctions. 
It can be seen in Fig. [I] that two trans and two cis min- 
ima exist in the calculation. Each well possesses a twofold 
symmetry, such that the PES has four equivalent (trian- 
gular) quadrants. More generally, the symmetry group of 
the reduced dimension Hamiltonian is of order four. Re- 
flection across the antidiagonal corresponds generally to 
the (12) (a6) CNPI symmetry operatior^, but can also be 
thought of as the operation in the C2V point group. 
Reflection across the diagonal is similarly (12)(a6)*, or 
the i operation in the C2h point group. This leads to 
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TABLE IV. Rovibrational symmetries of DVR wavefunctions 



B. Agreement with Experiment 



CNPI-MS^ 


C2h 


C2v 


Sa - /As- 


ag, 


odd 


62 00 


Ss + /Aa+ 


ag, 


Kc even 


Hi ee 


Ss - /Aa- 


bu, 


odd 


ai 00 


Sa + /As+ 


bu, 


Kc even 


62 ee 



is not an equivalent rotation for the trans geometry, and 
therefore its oe/ee and 00/eo have the same CNPI symmetry. 
Conversely, the cis rotational structure supports all CNPI 
symmetries, but only 00 and ee appear here. The first member 
of a pair of CNPI-MS labels will be used as a shorthand 
notation. 



the conclusion that inversion through the center of the 
coordinate space as depicted in Fig. [1] correlates with the 
E* operation. Therefore, the eigenstates produced by 
the calculation will belong to one of the four rovibra- 
tional irreducible representations laid out in Table ITVl as 
will be illustrated in Fig. [2] Accordingly, we expect each 
calculated vibrational level to appear as a near degener- 
ate pair of eigenstates with different rotational symme- 
tries, except when cis-trans interaction causes one mem- 
ber of the pair with a particular rovibrational symmetry 
to shift. This shift is analogous to the iiT-staggering that 
occurs in the experimental spectrum^ ^. 
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FIG. 2. Wavefunctions illustrating the 4 symmetries listed in Table 
llVl and the types of levels discussed in Sec. lIII Bl C. Clockwise from 
the top left: trans 3^, Sa—; trans 6*, Ss-i-; trans 3^6^, Ss—; cis 6^, 
Sa-i-. These states belong to the classes of levels with excitation 
in only totally symmetric modes, pure bending levels, stretch-bend 
combination levels, and cis levels, respectively. Note that the cis 
6^ wavefunction is delocalized over both cis and trans wells. The 
nodal patterns for trans 113 are similar to those obtained by spectral 
quantizationi^ii&. 



The primary difference between the DVR results and 
the observed levels is that many states are missing due to 
the reduced dimensionality. Most important, due to the 
exclusion of V4, all bending (B) polyads are represented 
by a single vibrational level, i.e. instead of the n+1 vibra- 
tional levels in a i?" polyad, the 3D calculation includes 
only one level, 6". States involving quanta in vi and U5 
are also missing, but less conspicuously. 

The simplest illustration of the agreement between cal- 
culation and experiment is the juxtaposition of one pho- 
ton spectra in Fig. [3l (The method for computing spec- 
tra is described in Appendix A.) The correspondence be- 
tween individual features is clear at low vibrational en- 
ergy, not only for the strong Franck-Condon active pro- 
gressions that can be seen easily in the figure, but also for 
the weaker bending polyad o^'^" that have been identified. 
The comparison does become slightly more complicated 
at higher energy where the experimental characteriza- 
tion of the level structure is less complete. The onset of 
predissociation^^ above 46074 cm"-'^ causes many states 
to essentially disappear from the LIF spectrum, but the 
simulated spectrum makes no allowance for this effect. 
Furthermore, the density of states not present in the 
reduced dimension calculation increases rapidly above 
the fundamentals of the neglected CH stretching modes 
(45077.65 and 45054.97 cm~^ for vi and W5, respectively). 
Fortunately, these discrepancies do not cause significant 
difficulties in analyzing the calculated levels and extract- 
ing information relevant to the states observed experi- 
mentally. 
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FIG. 3. Comparison of calculated spectrum (downward) with 
the one photon LIEiS. The more intense peaks have been trun- 
cated to make the weaker transitions visible. Members of the two 
main Pranck-Condon active progressions are identified for refer- 
ence. Three arrows mark the approximate positions of (from left 
to right) the cis origin, the C2H-I-H dissociation limit, and the zero 
point energy corrected cis-trans barrier. 

We now proceed with a more detailed comparison 
of the calculated level structure to that observed ex- 
perimentally. In the recent experimental literature, 
the vibrational states are categorized in three groups: 
states containing excitation only in the totally symmetric 
modes^^'"*^, the pure bending polyads^, and the stretch- 
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bend combination polyadsiS. We organize our discussion 
accordingly. 



1. Totally Symmetric Modes 

The states 2"3™ are relatively straightforward, and the 
calculated levels match well with experiment. Values ob- 
tained from fitting the groups 2" and 3" agree with Table 
5 of Rcf. 30], with ijJ2 and approximately 20 cm~^ too 
high and too low, respectively. The parameters 0:22 and 
a;33 are close to correct but too small in magnitude, such 
that the residual for changes sign at Zv^ . Even for V2 , 
the residuals remain in the tens of cm~^ over the energy 
region of interest. The deviations from experiment are 
presumably due primarily to the neglect of the symmet- 
ric CH stretch. 



2. Pure Bending States 

We next consider the pure bending polyads^, where the 
lowest members are nominally the 6" states that exist 
in the DVR calculation. Here the comparison is compli- 
cated by the omission of the torsion from the calculation, 
since and W4 interact very strongly via Coriolis effects 
and Darling-Dennison resonance. Nevertheless, we find 
that loq from the DVR is only 8 cm"-'^ higher than the 
experimental value, and the anharmonicity is once again 
too small in magnitude. This is undoubtedly due to the 
surprisingly large contribution of the vibrational angular 
momentum to the experimental value of xge, an effect ab- 
sent from the calculation. Overall the bending behavior 
here is nearly harmonic, as is the case with the deper- 
turbed level structure. As in the previous section, the 
too high harmonic frequency is presumably due at least 
in part to the neglect of the CH stretches, since is also 
of bu symmetry. 



3. Stretch-Bend Combination Polyads 

The stretch-bend combination polyads demand a more 
detailed comparison, because a global model that ac- 
counts for their vibrational structure has not yet been 
developed, despite the existence of extensive assignments 
and rotational analysea^S. Such a comparison is pre- 
sented graphically in Fig. U] by plotting the effective fre- 
quencies of the bending modes as a function of quanta 
in vj,. It can be seen that the calculation reproduces 
well even the more unusual features of the observed level 
structure, and the minor differences are due to the dis- 
agreement in the diagonal anharmonicities, explained in 
the previous two sections. The pathological behavior of 
this set of states is not entirely unexpected, as a combina- 
tion of D3 and essentially constitutes the isomerization 
path coordinate. Excitation in both these modes in ei- 
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FIG. 4. Plots oiuj'^ and lOf^ vs. 1)3 derived from the observed 
and calculated level structure. Different marker styles denote 
data from different values of ug. Deperturbed values for To are 
used, except for the 3''6^ and S^G'^ levels, which have not yet 
been deperturbed. (In those cases the energy of the observed 
J = if = state is used instead.) 



ther well should promote cis-trans tunneling, which will 
cause mixings as discussed in subsequent sections. 

One oddity of particular interest is the dramatic de- 
crease in ojg^ for the 3"6^ series, illustrated in Fig. 13 of 
Ref. [l3 up to 3^6^. Our recent high sensitivity spectra 
have revealed a band that is a promising candidate for 
3'^6^, although the assignment has yet to be confirmed 
by rotational analysis and therefore details regarding it 
will be communicated later. We nevertheless incorporate 
it in the figure for comparison with the DVR prediction. 
We further include the well-known level at 47206 cm^^ 
as 3'^6^, based on the preliminary assignment of 3^6^ as 
well as previous discussions^Sii^, in addition to the strong 
intensity and proximity of 3*6^ to 3^ in the calculated 
spectra. For both of these tentative assignments we find 
good agreement with the DVR results. 

The importance of this sharp decrease in the effective 
frequencies is that it signals the onset of the cis- trans iso- 
merization process as the potential softens approaching 
the transition state. We now turn our attention to the 
cis minimum of the 5*1 state and its vibrational levels. 
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C. cis states 

The adiabatic energy separation between the cis and 
trans isomers of acetylene is found to be 2664 cni~^ 
based on high-level ab initio calculations. Specifically, 
the effects of valence electron correlation have been in- 
cluded up to the coupled-cluster singles, doubles, triples 
and quadruples (CCSDTQ) level of theory, using the 
equation-of-motion (EOM) variant of coupled-cluster 
theory to treat these excited state isomers. In addition, 
effects of basis set insufhciency are estimated using ex- 
trapolation techniques^, and contributions due to zero 
point energies, core correlation and scalar relativistic ef- 
fects are included as well. The final result of these cal- 
culations predicts the zero-point level of the cis isomer 
to lie at 44861 cm~^, a calculation that we believe to be 
in error by no more than 50 cm~^. In passing, we note 
that the present results are in line with a similar estimate 
published some time ago by Kallay and Gauss^° of 44852 
cm~^. 

At the EOM-CCSDT/ANOl level of theory used in 
the DVR, the cis ground state lies above that of the 
trans conformer by '^2820 cm~^, and its vibrational state 
manifold is consequently less dense at any given energy. 
We find that all but the lowest-lying states contain at 
least a few percent trans character, using the crude met- 
ric 2 Jq J^^ IV'P dRccd4>id(j}2, but even these states 
are predicted to have non-negligible intensity in one of 
the four spectra discussed in Ref. These results are 
partially summarized in Table [V] At higher energy it is 
frequently the case that one or both rovibrational sym- 
metries interact so strongly with several trans states that 
it is difficult to attach that zero-order assignment to any 
particular eigenstate. Accordingly, it should be kept in 
mind that the mixing fractions are probably sensitively 
dependent on local resonances. Spectral intensity is, as 
expected, generally correlated with trans character, and, 
as noted earlier, excitation in V3 and vq greatly enhances 
cis-trans mixing. An in depth investigation of cis-trans 
interaction and intensity borrowing for a specific level is 
undertaken in the next section. 



D. Investigation of a Specific cis -f^ trans Interaction 

In order to determine the effects of cis states on the 
trans level structure and possible sources for their inten- 
sity, we would ideally like to compare the predictions of 
a model that neglects cis-trans mixing to the true spec- 
trum (either calculated or experimental). In the absence 
of a global H°*^ for the Si state, it is difficult to con- 
sider the level structure and spectral intensities in terms 
of a zero-order picture perturbed by the addition of an 
interaction. However, we can use the DVR to approxi- 
mate two non-interacting minima by performing two cal- 
culations wherein the wavefunctions are restricted to one 
geometry or the other. The cis-trans "interactions" are 
then calculated by the trick of using the eigenstates of 



TABLE V. Calculated as levels up to 5,000 cm"^ 
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^ Values given are averaged between 00 and ee. 
This level will be discussed in detail in Sec. IIII Dl 

the full PES, I'ipi), as a complete set to find the vibra- 
tional overlap integrals between the zero-order cis and 
trans wavefunctions 

i 

and the mixing fractions follow from dividing by the zero- 
order energy differences. (The "true" interaction matrix 
elements also include an electronic pre- factor, so that we 
obtain only relative coupling strengths.) Zero-order spec- 
tra can also be computed between the X state and the 
non-interacting sets of Si states, which allows for the 
tracking of intensity borrowing when the interaction is 
turned on. It should also be noted that this method can 
aid in assigning highly mixed eigenstates, particularly by 
inspection of the quantities {ipi\<fm^)^ and {^pi\(p'^°''^'')° . 

The results of the above procedure as they apply to cis 
3^6^, a state that lies in the region below trans 3*, are 
discussed in the remainder of this section^^. The vibra- 
tional 62 symmetry of this state means that its if = 1 
level interacts with K — 1 ag levels, and that its K — 0,2 
levels interact with K — 0,2 6„ levels. The two rotational 
symmetries therefore require parallel but separate anal- 
yses, and so we will treat only 00, since it appears in the 
simpler one photon spectrumii.. 

The zero-order state of interest is trivially recognized 
by its nodal pattern (Fig. [5^) . The eigenstate with this 
nominal assignment can then be identified (Fig.[SjD), even 
though the zero-order state is mixed into several eigen- 
states. The spectra for the zero-order states and the 
eigenstates are plotted in Fig. ^-d for 100 cm~^ above 
and below cis 3^6^. These spectra differ in two diag- 
nostically important ways. First, the intensity has been 
redistributed such that the cis state has increased its in- 
tensity from essentially nothing to greater than that of 
trans 2^3^6^. This tells us not only that there is cis-trans 
interaction, but also that the intensity of the cis state 
derives entirely from this mixing and not from purely 
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TABLE VI. Ab mitio geometry and harmonic frequencies for 
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The transition state geometry is sliglitly cis-bent. 



vibronic effects, such as the variation of the electronic 
transition moment with the nuclear coordinates (Fig. [7]). 
Second, the energy level pattern has changed because 
the interacting states repel one another. Interestingly, 
although ordinarily no strong resonance would mix trans 
2^3^6^ and trans 3*, their interactions with a common cis 
state cause them to move apart, an example of indirect 
mixing. 

The explanation for these two phenomena is displayed 
in Fig. [5]3-f. We see that trans 2^3^6^ has a much larger 
overlap integral with cis 3^6^, as expected because of its 
excitation in wg, but trans 3* is near resonant, and conse- 
quently the two states have approximately equal mixing 
angles with cis 3^6^. This provides the cis state with 
significant intensity, and in fact the calculated relative 
intensities for these three states agree quite well with 
those observed experimentally. Finally, the vibrational 
overlaps between the two trans states and the cis state 
are in fact of opposite sign, which suggests the possibility 
of observable interference effects. This three state inter- 
action is currently being investigated experimentally, and 
early results are so far consistent with the calculation. 



E. Above-Barrier Dynamics 

A full dimensional transition state search finds the cis- 
trans barrier height to be 4979 cm~^ (Table IVip . In 
the reduced dimension PES used in the DVR, the bar- 
rier height is effectively 5145 cm~^i^ Above this energy, 
delocahzed states begin to appear— (Fig. but other 
states unrelated to the isomerization coordinate remain 
unaffected. Although any conclusions about the above- 
barrier dynamics would be premature, a cursory analysis 
of the delocahzed wavefunctions shows that they corre- 
spond to out of phase combinations of local benders. The 
elliptical shapes ( "ring modes" ) are therefore similar to 
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FIG. 5. Plots depicting the analysis of cis-trans mixing 
and intensity borrowing involving cis 3^6^ oo. (a) Zero-order 
wavefunction of cis 3^6^ oo. (b) Eigenstate with nominal 
assignment of cis 3^6^ oo, showing significant trans charac- 
ter, (c) Stick spectrum of the zero-order states. The peak 
of cis 3^6^ 00 (white) has been multiplied by 10^. (d) Stick 
spectrum of the eigenstates, in which intensity has been trans- 
ferred to cis 3^6^ oo from the nearby trans states. (The inten- 
sity of trans 3* is off scale in both (c) and (d), and the energy 
offset between the two spectra is slightly arbitrary.) (e) Mag- 
nitudes of vibrational overlaps between the states shown in 
(c) and cis 3^6^ oo. Many other states (not shown) have as 
large or larger vibrational overlaps with cis 3^6^. (f) Mag- 
nitudes of the mixing angles between the states shown in (c) 
and cis 3^6^ oo. These are the largest mixing angles by at 
least an order of magnitude, belonging to trans 2^3^6^ and 
trans 3^, whose zero-order wavefunctions are shown in (g) and 
(h), respectively. 



Lissajous figures for two equal amplitude oscillations with 
a phase difference of ±7r/4 or ±37r/4, unlike the patterns 
along diagonal or cardinal axes that accompany normal 
or local mode behavior, respectively (Fig. The two 
differently inclined ring modes are presumably the above- 
barrier counterparts of cis and trans bending. 

Some states that have amplitude in both wells do not 
exhibit this ring mode behavior; in our results such states 



9 





150 200 
ZCCH1 



1 50 200 250 
iCCHI 



the calculated results to investigate cis-trans mixings. 
Although we believe that the reduced dimension DVR 
method demonstrated here is a powerful approach for 
studying isomerizing systems, there are interesting as- 
pects of the Si C2H2 spectrum that it does not address, 
including strong vibration-rotation interactions and pos- 
sible multiple pathways to isomerization. We intend to 
develop a full dimensional treatment for Si acetylene to 
study these effects and enable comparison with the com- 
plete set of experimental observations. 



FIG. 6. Wavefunctions of delocalized states above the iso- 
merization barrier. 

always have nodal patterns in each well that are clearly 
assignable, indicating that these states are less affected 
by the isomerization despite their delocalization. 

The shapes of the wavefunctions both above and below 
the barrier to isomerization provide some clues about the 
important resonances in the bending Hamiltonian. Since 
we do not observe the cross shaped wavefunctions indica- 
tive of local modes, we can infer that the i^aaee Darling- 
Dennison resonance does not dominate the dynamics. 
This is in line with the frequency ratio of the two modes 
being not 2 : 2 but rather ~ 1.3 : 2, lying in between 
the usual ratios for strong Darling-Dennison and Fermi 
resonance. We find^ that both types of resonance are 
necessary in an H''*^ to reproduce the unusual nodal pat- 
terns of the below-barrier trans well DVR wavefunctions. 
The presence of both strong Fermi and Darling-Dennison 
resonance leads to the destruction of any polyad struc- 
ture in the A state. For fixed values of (wi, D2, all 
states of a given symmetry can interact via known anhar- 
monic interactions, so the Hamiltonian does not readily 
block diagonalize according to conserved polyad quan- 
tum numbers. Research into this phenomenon and its 
relation to low barrier isomerization is ongoing. 



IV. CONCLUSION 

We set out in this paper to investigate the spectro- 
scopic consequences of low barrier cis-trans isomerization 
in Si acetylene by calculating the vibrational eigenstates 
of a high accuracy PES using a reduced dimension DVR 
method. The calculation reproduces some difficult as- 
pects of the trans conformer level structure, which are ul- 
timately due to the isomerization. Another consequence 
of the isomerization process is that nominally forbidden 
transitions to the cis conformer appear. These transi- 
tions occur in our calculation near the observed energies 
of the "extra" levels in the A ^ X spectrum of acety- 
lene. We find that our reduced dimension DVR calcu- 
lation agrees with the explanation of these "extra" lev- 
els as belonging to the cis conformer, and that we are 
able to explain the intensities of these levels by using 
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Appendix A: Computation of Si o So Spectra 

In order to calculate A ^ X spectra, two additional 
quantities are required beyond the A state vibrational 
wavefunctions: the electronic transition moment, ^, and 
the X state vibrational wavefunctions. We calculate the 
transition moment ab initio over our coordinate grid at 
the EOM-CCSD/ANOl level of theory, again using the 
CFOUR program system''^. The X state vibrational 
wavefunctions can be obtained merely by replacing the 
A state PES with that of the X state and repeating the 
DVR calculation. This portability of the method is one of 
its powerful features. Although the results thus obtained 
for the X state are not worth examining in depth, given 
the numerous ab initio treatments in the literature, we 
would like to note that the qualitatively important fea- 
tures related to large amplitude motions are simulated 
correctly by our DVR method. Chief among these is the 
emergence of the "local-bender" states^ii^ from the nor- 
mal modes and v'^, as shown in Fig. [5] 

With these two results in hand, we calculate for 
each pair of eigenstates the spectral intensity /„m — 
|(m|/i|n)p, where m and n belong to different electronic 
states. In order to calculate the product it is necessary 
to interpolate the transition moment and wavefunctions 
to a common set of grid points. This process yields 
all possible upward spectra, including the four necessary 
to sample experimentally the full rovibrational structure 
of the A state^, and also all possible downward spectra 
(DF/SEP). This wealth of information was of enormous 
value when analyzing the calculated results, and we ex- 
pect it to guide future experiments. 
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FIG. 7. The transition moment /i as a function of the bend- 
ing angles. Note that /x = at all points of C2v symmetry, in- 
cluding linearity. In order for the overall symmetry of to be 
Ss—, the quantity used to calculate yl •«->■ X spectra must have 
Sa— symmetry in the DVR coordinate space. Briefly, this re- 
quirement arises because the X state DVR wavefunctions are 
vibrational instead of rovibrational. Adding the rotational 
symmetry factor necessary to ensure that transitions occur 
between the same nuclear spin symmetry species includes an 
extra Sa— in /i. 
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FIG. 8. From top to bottom, pairs of wavefunctions de- 
picting the evolution from the normal mode limit to the lo- 
cal mode limit in So acetylene from approximately 8, 000 to 
12,000 cm~^ of vibrational excitation. 
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